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Abstract 

Generalised Hagedorn wave packets appear as exact solutions of Schrodinger 
equations with quadratic, possibly complex, potential, and are given by a poly¬ 
nomial times a Gaussian. We show that the Wigner transform of generalised 
Hagedorn wave packets is a wave packet of the same type in phase space. The 
proofs build on a parametrisation via Lagrangian frames and a detailed analysis 
of the polynomial prefactors, including a novel Laguerre connection. Our find¬ 
ings directly imply the recently found tensor product structure of the Wigner 
transform of Hagedorn wave packets. 
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1. Introduction 

Hagedorn wave packets are prototypes of quantum wave functions that are 
highly localised in both, position and momentum space. They arrive as eigen¬ 
states of general multidimensional harmonic oscillators. The fc-th wave packet, 
k £ N d , can be written as 

<Pk{x) = Pk(x)g(x ), x £ R d , 

where pk is a polynomial of degree |fc| and g is a Gaussian with complex width 
matrix; see flOl fl2j [18j . 
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The favourable properties of Hagedorn wave packets have turned them into 
an essential tool for various applications, notably in molecular quantum dy¬ 
namics and quantum optics. For instance, they have been used in the analysis 
of time-dependent semiclassical Schrodinger equations nu m ns, the con¬ 
struction of numerical methods for high frequency problems DUE], and the 
representation of quasiprobabilities in quantum optics Ell¬ 
in the literature, Hagedorn wave packets also coexist as generalised squeezed 
states; see (5j. These states are generated by applying translation and squeezing 
operators to the eigenstates of the standard harmonic oscillator. However, this 
construction and Hagedorn’s approach produce, up to a phase factor, the same 
wave packets, see |[IB] . 

In this paper our aim is to analyse the representation of Hagedorn wave pack¬ 
ets on phase space, that is, their Wigner functions, and give a comprehensive 
description of the appearing structures. 

So far, Wigner functions have been classified for Hermite functions, see [7j 
120j , by using the well-known Hermite-Laguerre connection. This link has been 
generalised in m , showing that the Wigner functions of a Hagedorn wave packet 
factorises regardless of the structure of the wave packet itself. Another approach 
using ladder operators in the two-dimensional setting can be found in |.4j. 

Building on the generating function of the wave packets from m , we are 
able to present a general formula for the Wigner functions of Hagedorn wave 
packets. As our main result, in Theorem ED we show that the wave packets in 
phase space are given by Hagedorn wave packets of doubled dimension. This 
insight allows to lift all qualities of the wave packets to phase space and, thus, 
enables new perspectives and approaches in modelling and approximation the¬ 
ory. We find that the structure of the Wigner functions can be explained with 
the properties of the wave packets. A connection between the polynomial pref¬ 
actor in the Wigner functions and Laguerre polynomials then reproduces the 
previously found factorisation results. 

Wigner functions play a central role in microlocal analysis and semiclassical 
quantum mechanics; see, e.g., EMU- In particular, various important semiclas¬ 
sical approximations and algorithms for the simulation of quantum molecular 
dynamics combine the Wigner function of the initial state with classical dynam¬ 
ics on phase space, see, e.g., mmm- Our results can contribute to improv¬ 
ing these methods by providing new methods for the approximation of Wigner 
functions of more general states, and extending the applicability of wave packet 
propagation methods to phase space. 

1.1. Outline 

Motivated by the form the wave packets attain under non-unitary time evo¬ 
lution, in Chapter 2 we start with a generalised definition for Hagedorn wave 
packets; see also HU- For our analysis we follow a geometric approach based 
on Lagrangian frames. However, all conditions can easily be translated into the 
classical ones found in m or Ga¬ 
in order to find criterions for tensor factorisation, in Chapter 3 we analyse 
the polynomial prefactors pk of the wave packets. We characterise these poly- 
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nomials by a three-term recursion relation, a generating function as well as a 
ladder operator. All definitions are closely related to the corresponding formu¬ 
las for Hermite polynomials, and thus suggest an interpretation as multivariate 
Hermite polynomials. However, the polynomials are not simple tensor prod¬ 
ucts of univariate Hermite polynomials, but generically exhibit a more complex 
structure, as shown in Proposition [8] Nevertheless, our findings can be read as 
reasonable multivariate extension of the relations between Hermite and Laguerre 
polynomials discussed in [2UJ]. 

In Chapter 4, we use the generating function to identify the Wigner func¬ 
tion of Hagedorn wave packets with a wave packet of doubled dimension. This 
result combined with our polynomial analysis from Chapter 3 give a detailed 
picture of Hagedorn wave packets in phase space that constitutes the core of 
our manuscript. 

Finally, in Chapter 5, we illustrate characteristic examples for both, polyno¬ 
mials and wave packets in two dimensions. 

2. Generalised Hagedorn wave packets 

In this chapter we adopt the viewpoint of m, and briefly sketch the geo¬ 
metric approach to defining Hagedorn wave packets. Moreover, we extend this 
framework by decoupling the raising operator from the ground state. This type 
of generalised wave packets appears for example in the context of non-selfadjoint 
evolution problems, see 

2.1. Lagrangian frames and ground states 

We consider the classical phase space T*M. d = R 2d , equipped with the stan¬ 
dard symplectic form 

Definition 1. A matrix Z £ C 2dxd is called a Lagrangian frame, if 

z T nz = o (i) 

and rank Z = d. Furthermore, Z is called normalised, if 

l -z*nz = id d . ( 2 ) 

In the following, we will denote typical phase space points by z = ( q,p ) £ 
K 2<i , with q £ M. d being the position, and p £ M. d the momentum variable. 
Similarly, we write Z = ( Q;P ) with Q,P £ C dxd for Lagrangian frames. With 
this denotation, the isotropy condition ([l]) and the normalisation condition Q 
coincide with the symplecticity condition in m or, by taking Q = A and 
P = iB , with Hagedorn’s original definition in [T2| . 
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The image of a normalised Lagrangian frame Z is a complex Lagrangian 
subspace range Z C C 2d . All Lagrangian frames spanning the same subspace L 
are related by unitary transformations, see, e.g., m § 2 . 1 ], and thus have the 
same Hermitian square ZZ*. Hence, considering the matrix 


G z = n T Re(ZZ*)fl 


( PP* - PQ* + ild d \ 

\-QP* - ild d QQ* ) 


( 3 ) 


is a convenient way to classify the Lagrangian subspace range Z. The matrix 
Gz is real, symmetric, positive definite and symplectic, and thus called the 
symplectic metric associated with the Lagrangian subspace range Z. The matrix 
Gz plays an important role when lifting the wave packets to phase space. 

Following the idea of Hagedorn, we can associate a coherent ground state 
with any normalised Lagrangian frame. 

Lemma 2 (Ground state). Let Z = (Q;P) be a normalised Lagrangian frame. 
Then, Q and P are invertible and 

Im (PQ _1 ) =(QQ*)~ 1 >0. 

In particular, for e > 0, 

<Pq[Z](x) = (7re)“4 det(Q )“2 exp (^x T PQ~ 1 x) (4) 

is a square integrable function with ||[^]||^2 = 1. 

A rigorous proof of this result can for example be found in m ■ Briefly, 
isotropy ([l]) ensures <p%[Z\ £ L 2 (R d ), while (Pq[Z] is normalised if and only if Z 
is normalised. 


2.2. Excited states and spectral properties 

Let q denote the position operator and p the momentum operator. Analo¬ 
gously to phase space points, we write S' = ( q,p ). Using this notation, we can 
define a linear operator 

Al l Y } = 7Te Y * n * ( 5 ) 

associated with a normalised Lagrangian frame Y. The components of [U] 
commute due to the isotropy of Y . Starting from a ground state [Z\, the 
generalised Hagedorn wave packets ip e k [Z, Y] then are constructed via \Y\ as 
follows. 

Definition 3 (Generalised wave packets). Let k £ N d and Z,Y £ C 2dxd be 
normalised Lagrangian frames. Then, the k-th Hagedorn wave packet is defined 
as 

rt[Z,Y\ = j a (AW<fi[Z\ ( 6 ) 

with standard multiindex notation, (A^\Y]) k = AlfF ]* 1 ... A^\Y] k d d . 
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Based on this definition we will refer to [Y] as raising operator. 

If Y = Z, the above definition yields the standard Hagedorn wave packets 
{<p|[Z]} fcgN d that form an orthonormal basis of L 2 (M. d ). In this case, the adjoint 
operator of A'[Z\, 

MZM[Z) = Vk^i_ e .[z}, 

plays the role of a lowering operator. This property does not hold for the 
generalised wave packets. 

The construction of Hagedorn wave packets by means of A^[Y] implies 


<Pl[Z,Y] = pF k [Z,Y]<p* 0 [Z), 

where p\[Z,Y], k £ N d , is a polynomial of total degree \k\. 

Proposition 4. Let Z = (Q;P) £ C 2dxd and Y = (X;K) £ C 2dxd be two 
normalised Lagrangian frames and B = ^Z*YIY. Then, for k £ it holds 

Vl{Z,Y}(x) = (^YB*Q~ 1 x s J Vo[Z](x) (7) 

with M = jY*GzY + B*Q~ X QB and the polynomials {q k d jfegN" 1 that are re¬ 
cursively defined by q q 1 = 1 and 

(Qk + e 3 (x))U = - 2M ■ (kjQk-^U- 


In the special case Y = Z we find B = Id and M = Q l Q. 

Proposition [4] is easily proven by identifying the found raising operator with 


the recursive definition of {v 3 *;} fc? see Appendix Appendix A A similar result for 


special choices of Y can be found in [17. . Proposition|4|reveals that the structure 
of the Hagedorn wave packets can be understood by studying the polynomial 
prefactors {q k r }k- This polynomial analysis is conducted in the next section. 

In this chapter we only introduced wave packets centered at the origin. In 
general, one obtains wave packets centered at any phase space point z = (q,p) £ 
R 2d by applying the linear Heisenberg-Weyl operator 

(T z i/>)(x) = e ipT{x - q/2) /^(x - q) 


to the corresponding wave packet centered at the origin. 

Also in the remaining parts of this paper, we only present our analysis for the 
wave packets centered at the origin. However, we stress that all of our results 
and proofs easily adapt to the more general case. 


3. Analysis of the polynomial prefactor 

In this chapter we analyse the polynomial prefactors of Hagedorn’s wave 
packets, which are defined for a symmetric and unitary matrix M £ C dxd via 
the three-term recursion relation (TTRR) 

(9& e >))i=1 = 2x ^ (*) - 2M ■ (Mfe»);=l, (8) 
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see Proposition [ 4 ] where denotes the j-th unit vector in with boundary 
conditions q q 1 = 1 and qf 1 = 0 for all l ^ This recursion is well-defined as 
M is symmetric. Equivalently, the polynomials q^ 1 could be defined via their 
generating function or raising operator, both of which we derive in this section. 

Surprisingly, this class of polynomials seems to be little studied in the lit¬ 
erature found by the authors. A noteable exception is [5] who defined them 
through their generating function and showed a generalised Mehler formula. In 
jl] one can find results for the two-dimensional case. 


3.1. Generating function and raising operator 

To get an intuition, we quickly discuss the univariate polynomials first. In 
one dimension ([8]) simplifies to 

H*+ 1 (®) = 2 xH*(x) - 2XnH^_ 1 (x) 

with A £ M. Starting from Hq = 1 and Hf 1 = 0, for A = 1 we produce the 
usual Hermite polynomials. For A ^ 0 we find rescaled Hermite polynomials 


Ht = \'2 Hi 



To complete the picture, for A = 0, we generate monomials H® = ( 2x) n . 

In the multivariate case the polynomials may emerge as simple tensor prod¬ 
ucts of Hermite polynomials, but typically attain a more involved structure. 
Recently, in m, there has been published a formula for the generating function 
of (rescaled) polynomials of type (j8j. 

Lemma 5 (Generating function). Let M £ <C dxd be symmetric and unitary. 
Then, the generating function of the polynomials {Q^}keN d given by 


j.k 

f(x,t)= ^2 ff q k(x) = exp(2x T t - t T Mt). (9) 

fee N d 


To give a more profound characterisation of the polynomials, we also look 
at the raising operator and the gradient formula for the polynomials q^ 1 . In the 
language of Dirac ladders, the gradient plays the role of a lowering operator or 
annihilator for the polynomials. 

Lemma 6 (Ladder operators). Let b' M = 2x — MX7 x and k £ N d . Then, it 
holds 

(q^+e j )U= b Uf and V 9 f = 2(k jq ™ e .)U- ( 10 ) 

Proof. The generating function satisfies X7 x f = 2tf and 

V*/ = 2 xf - MV x f = tf M f- 

By the definition of the generating function f(x,t) it follows that 


d t] f(x,t) = 

k 


t k ~ e 3 

{k-ej)\ 
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i.e. dtjf simply shifts the summation index by one into the direction of ej. The 
gradient formula then follows immediately from (|8|. □ 

As a consequence of (|8j) and Lemma[6] the polynomials satisfy the Rodrigues 
formula 


qjf(x) = exp (x t M~ 1 x) (-AfV) fe exp (-x T M“ 1 x) . (11) 


3.2. Factorisation and eigenvectors 

The structure of the multivariate polynomials qjjf crucially depends on the 
matrix M in the TTRR ([8]). In fact, the polynomials qjf factorise for all k £ N d 
if and only if their generating function factorises. We summarize this observation 
in the following remark. 

Remark 7. By inspecting Lemma [5] one observes that the generating function 
of the polynomials factorises into m lower-dimensional generating 

functions of the same type if and only if there is a relabeling of the coordinates 
such that M is block diagonal with m blocks. In this case, the polynomials 
qjY are tensor products of m lower-variate polynomials of the same type for all 
k £ N d . 

In particular, the polynomials qj are tensor products of univariate polyno¬ 
mials for all k if M is diagonal. On the contrary, if M has nonzero offdiagonal 
entries, the polynomials q£ r are not simple tensor products anymore. 

The polynomials qjY form a set of simultaneous eigenvectors for the operators 


Tj = 


(6 f 


MJJ 


)jd Xj + d x , {b\j)j 


= (1 


2 Xjd x 


)-cL.(MV),- , j = 1.rf. 


acting, for example, on the space of polynomials on R d . This can easily be seen 
from 


T n M — I 
1 J < ik 2 


\.Ok n M 


a M 

jQk+e 


— (2 kj + 1 )q 1 jf , 


similarly as in the computation of harmonic oscillator eigenvalues. Hence, qjf 
is an eigenvector belonging to the eigenvalue 2 kj + 1 of Tj. Moreover, for j ^ k, 
the operators Tj and T ^ commute if and only if Mjk = = 0. In other words, 

by the remark above, commutation of the operators reflects the factorisation of 
the polynomials. 


3.3. Laguerre connection 

As noted in Remark [ 7 ] the polynomials qjf factorise for all k if and only if 
M is block-diagonal. Therefore, our aim is to express the general polynomials 
as a linear combination of tensor products by deleting off-diagonal entries of M. 
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It turns out, that q™ can be rewritten as a Laguerre polynomial 


4 a) (*) = E 

3= 0 


— x) ] , n £ N , a > 0, 

n-j ) j\ y 


of the raising operators corresponding to the polynomials generated by a reduced 
matrix applied to 1. This generalised Laguerre connection adds a new point 
to the long list of relations between Hermite and Laguerre polynomials, see, 
e.g., M- More precisely, we express qff via raising operators where 

M[n, to] denotes the matrix M with deleted offdiagonal entries M nrn and M mn , 
i.e. 

0 , = , 


(M[n, m]) ■ ■ = 


Mij , otherwise. 


Proposition 8 (Laguerre connection). Let M £ C dxd be symmetric, and 
M nm = A jtz 0 for some n ^ m. Suppose k £ N d with k n > k m . Then, 


Qk 


(*) = (4*) 


l = U\ k - k m{e n +e m ) 


{-2\) k -k m \L^~ km) (£440 i, 


where c* = b' M ^ n m j denotes the polynomial raising operator for the reduced ma¬ 
trix M[n,m\. The case k n < k m is analogous. 

Proof. Let /m denote the generating function of the polynomials qjf derived 
in Lemma [5] and define k[n,m] = k — e n k n — e m k m . Then, one computes the 
series expansion in t as 

f M {x, t ) =f M [n,m] (x, t) exp(—2Af„f m ) 

= ( E f7\ ( cf ) k ] ’ (! - 2A t n t m + jr(2A t n t m ) 2 - .■.) 


<ke N d 


= £ 




fee N d 


( k[n , to ])! 


E 

3=0 


j ! ( k n -j)\ {k m -j)\ 


Due to the definition of the generating function, this implies 
min(fe„,fe m ) If—2AV 

i: ’ ..n (4) fc -- j (4 t )*"- j 1 


^ rx 

,;'(x)=(ct)'= i " mi ^ 


J=0 


-j)K k m ~jV- 


( 12 ) 


and we can reorder the sum by means of the index t = k m —j > 0, since k n > k„ 
holds by assumption. This finally leads to 


{x) = (J) kln ’ m] (- 2 A) fc "*fe m !(4) fc "- fc ™ E 


e=o 


^n-\ 2\ n mJ _ 

{k m ~ P)Kk n ~ k m + ^)!^! 


1 


= ( c t) fc[n ’ ml (-2A ) k ~k m \{ci) k ~- k -L^ : - k ^ (^44.) 1, 

where we utilised that cL and cL commute. 


□ 








A direct, illustrative consequence can be deduced for the two-dimensional 
case, for which 


M = 


X\ A3 
A3 A2 


with A 3 , A 2 , A 3 £ C. 


Corollary 9. Let k £ N 2 with k\ > Then, 

'(—2A 3 ) fc 2 fc 2 !(ai i ) fel - fe 2 L^- fe2) (^ajj 1, A 3 7 ^ 0, 

K) fcl (4 2 ) fe2 i, 


= 


As = 0. 


(13) 


The case k\ < k 2 is analogous. 


For the case A 3 = 0, Equation (131 directly carries out to the factorisation 


In the case A 3 7 ^ 0, (131 guarantees that each qjf is just a linear combination of 
at most minjfci, A 2 } many tensor products of the form 

4J (4J 1 = H*(xi)H%(x 2 ) (14) 


where n — m = k\ — ^2 and k\ — k 2 < n < k\, m < k- 2 - Moreover, if Ai = 
A 2 = 0 , the diagonal creation operators aj produce monomials, and we obtain 
the formula 


<£{*) = {-\^k^x k r^L^\l 3 x lX2 ) (15) 

whenever A 3 7 ^ 0 and k\ > /c 2 - See also §5.1| for illustrations. 

We note that by applying Proposition [ 8 ] iteratively one obtains an expansion 
of the general polynomials in terms of tensor product Hermite polynomials, 
see Appendix |Appendix B| 


4. Hagedorn wave packets in phase space 

There are various representations of quantum systems on the classical phase 
space R 2d . The most popular one is the Weyl correspondence, where suitable 
selfadjoint quantum observables A : L 2 (M. d ) —7 L 2 (M. d ) are represented by their 
semiclassical Weyl symbol a : R 2d —> R, and wave functions £ L 2 (R d ) by 
their Wigner function 

= (2Tre)~ d f Jp(x + %)ip(x - |)e z? y/e dy. (16) 
J B d 

Then, matrix elements of A can be computed via the phase space integral 

(‘P, A tl>) L 2(Rd) = f W e (<p,ip)(z)a{z)dz. (17) 

J R2d 
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In this chapter our aim is to analyse the Wigner functions of Hagedorn wave 
packets. For readability, we write 


m,AZ, Y \ = W s mZ,Y\,vl[Z,Y}) 


for k,£ £ N d , and regard (k, t) as a multiindex in N 2d . Note, that by invoking [3] 
Equation (9.25)] we could also allow two wave packets that have different phase 
space centers. 

As our main result, we prove that the Wigner functions of Hagedorn wave 
packets on R d are given by related Hagedorn wave packets on the phase space 
R 2d . For proving this invariance result, we first introduce a phase space lift of 
Lagrangian frames. 


Lemma 10. Let Z = ( Q;P ) £ C 2dxd an d y = (X\K) £ C 2dxd & e i wo nor¬ 
malised, Lagrangian frames and B = | Z*QY. We define the lifted matrices 
as 


Z = 




and y 




in C 4dx2d with the symplectic form f l^ d = Q <g> Ido on R 4d . These have the 
following properties 

1. Z and y are normalised Lagrangian frames if and only if Z and Y are. 

2 . The symplectic metric fulfills VQY 1 = 2 iGz- 

3. For the lifted frames we have 


B = , z *n 4d y=( B Q ° B y 

4. The lifted recursion matrix satisfies 

(-\Y T G Z Y \ 

\ B*B lY T G z Yj' 

5. For the special case Y = Z it holds B = Id 2 d and 


M = Q~ l Q 


( 0 Id d \ 
\Idd 0 ) ■ 


Proof. For the first assertion, one computes 


y T tu d y 


(—Y T flY 0 \ 

^ 0 Y T flYj ’ 


y*n 4d y 


f-Y*flY 0 \ 

^ o y*qy)' 


and notes that y is normalised if and only if Y is. 

One can easily prove that Q -1 = iP*. Hence, the second part follows from 


VV = (S1Z - SI Z) ( /.(f) 


-n ( ZZ T + ZZ*) H = 2 Gz. 
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The formula for B is a direct computation. For the claimed form of M first 
compute 


Q~ X Q 


/- z T nz 

-z*nz\ 


ld d \ 

1 z*nz 

z T nz ) 

Udd 

o J 


By calculating 

y T c z y = y T nJ d zz*n 4d y 


f-Y T G z Y 0 \ 

'v 0 Y t G z y) 


the claim follows. The last case follows from the isotropy condition under the 
additional assumption. □ 


Since Z is again a Lagrangian frame, we can lift all our previous results to 
the phase space and consequentially find a family of Hagedorn wave packets 
in doubled dimension. In order to avoid confusion, we denote Hagedorn wave 
packets on phase space by upper case letters hr particular, a direct 

computation shows that 


W% 0 [Z])(*) = (• ne )- d e-* TG = (2tt e )- d / 2 S \ m [ Z ]{ z ), (18) 

see, e.g., m- For excited wave packets the following result holds true. 

Theorem 11. Assume that Z,Y £ C 2dxd are normalised Lagrangian frames. 
Then, for k, l £ the Wigner function Wk,e[Z,Y] is a Hagedorn wave packet 
on phase space, 

yvi l [z,Y] = {2ne)- d ^i k/) [z,y\. 

Consequently, it holds 


MA z ’ y ]( z ) 


(27re) d/2 M 
V2\ k \+\*\k\e'. e) 




where the lifted matrices y,Z , Q,V,B, and Ai have been defined in Lemma [7d| 


Proof. For the generalised Hagedorn wave packets, we find the generating func¬ 
tion 


E 

fceN d 


t k 

71 


72l7<4 [Z,Y](x) 


- 7 = t T B* Q ^X — t T Mt £\ry~\( \ 

e^ Vol z \( x ) 


h t {x) 


(19) 


by identifying the polynomial factors in Proposition [4] Note that h t is a Gaus¬ 
sian function in x and t. By writing z = (a;,£) and v = ( t,s ), we can easily 
compute the Wigner transformation of ht and h s as 


W £ (/i t ,/i s )(z) 


(2 tt£-) d [ h t (x+ ^)h s {x - dy 

Js. d 2 z 

(TT £r d e -z T G Z z/e e ^ vT B*Q-G+ vTMv ^ 
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where we identified the lifted matrices from Lemma [TOl 

Then, by formally interchanging the order of integration and summation, it 
is clear that 

W^h u h s ){z)= £ -1^21*1+1*1 W M [Z,y](z), 

MeN d 


which implies the result. For a rigorous justification, note that due to the 
Gaussian decay we can differentiate under the integral sign with respect to 
v = (t,s). Evaluating the differentiated function at v = (t, s) = (0,0) then 
shows the result. □ 


Remark 12. By invoking part (5) of Lemma |10| Theorem 11 explains the fac¬ 
torisation result for Hagedorn wave packets that has recently been discovered 
in HD. Namely, the polynomial prefactor of Wf. ( [Z, Z] is a product of d poly¬ 
nomials, 


mA z ’ z \( z ) 


(27 T£)~ d / 2 

V2\ k \+Wk\£\ 


*§[*](*> 


1=1 



0 V 

i or 


where, for kj > tj and N = 


Ut> lo T ( fe i ^ ) {2x 1 X 2 ) 


is a Laguerre polynomial of the form (15). 


5. Examples 

5.1. Polynomial prefactor 

In order to illustrate different types of polynomials qjf 1 , we present various 
examples for the nodal sets of two-dimensional polynomials qfj, for (i. j) £ N 2 . 
For simplicity, we restrict ourselves to real matrices M such that the polynomials 
generated by the TTRR (|8| have real coefficients. 

As examples we consider the unitary, symmetric matrices 

m<,, =G ;)• m<2, =(i ;)■ m<3> =*(; -O- <2o) 


By recalling and @, the polynomial q^g 1 corresponds to a simple 
tensor product of one dimensional Hermite polynomials, while 


= 6!2 7 x 1 L ( 6 1) (2x 1 x 2 ). 


( 21 ) 


One can see the consequences of these simple formulas in the structure of the 
nodal sets depicted in the upper panels of figure [T] 
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The matrix M^ gives rise to a more complicated mixing between the two 
variables. This can also be seen from the illustration of in the lower panel 
of figure |T| It is striking that already for real matrices M in two dimensions the 
polynomials generated by the simple TTRR ^ develop such nontrivial nodal 
sets. 



-5 0 5 




Figure 1: The nodal sets of three examplary two-dimensional polynomials, associated with 
the matrices Aft) (upper left), A/) 2 ) (upper right) and Aft 3 ) (lower). Regions with negative 
values are highlighted by grey coloring. 


5.2. Hagedom wave packets 

In this section we present two-dimensional examples of Hagedorn wave pack¬ 
ets in order to indicate the variety of structures that can be realised. 
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For our illustrations we employ the same matrices j £ {1,2,3} from (201 
as used for the polynomials in §5.1| Our choice for the Lagrangian frames 
Zj — ( Qj',Pj) satisfying 


Qj'Q, = M u \ j = 1,2,3, 

is given by 


1 

f 1 M 

1 -1 

„ 1 

/l + i 

1 — i 

i-zX 
1 + i 


( i 

1 

—i(l + v^2)\ 

y/ 2-1 

75 

i i 

’ 2 

i -1 

2 

i +1 

2 > 

, Z 3 = 

l—y /2 

2 \h 

1 

2 V 2 


v -v 


1+1 
\ 2 

w 1 ) 


2 

V 2 sh 

-J— 

2 V 2 / 


One can easily check that Zi,Z 2 , and Z 3 are normalised Lagrangian frames. 

We also consider an example of a generalised wave packet associated with 
the two Lagrangian frames Z 2 and Z%. In this case, the mixing matrix of the 
polynomials is given by the formula in Proposition |4| and does not equal Q~ X Q. 
We stress that despite the fact that the polynomials 

-Mj . 1 O Q 

Q k > J = 1,2,3, 

have real coefficients, the wave packets itself are not real-valued since the poly¬ 
nomials are evaluated on the subspace C C d , which does not coincide 

with except for Z\. This case is very special, since one has 

M = Id <*=*► Q £ R dxd , 


which implies that the standard tensor Hermite polynomials appear only 
together with real transformations Q. Hence, all Hagedorn wave packets with 
Q~ X Q = Id correspond to rescaled, sheared, or shifted multivariate Hermite 
functions, while this is not true in the case Q (j M. dxd . 

Figure[2]displays the absolute value of two Hagedorn wave packets associated 
with the Lagrangian frames Z\ and Z 2 . One can recognise that the Hagedorn 
wave packet <^\[Z\\ on the left hand side is just a rotated and rescaled Hermite 
function, as expected. In contrast, the wave packet l [Z 2 \ associated with 
Mt 2 ' 1 has a circular structure, which arises due to a complex rotation of the 
hyperbolas from the upper right panel of figure [T] 

For Z 3 the resulting wave packets exhibit complicated structures. The same 
is true for generalized wave packets associated with the two different Lagrangian 
frames Z 2 and Z 3 , as illustrated by the examples in figure [3] 

For us, the variety of different Hagedorn wave packets is very fascinating. We 
suggest that the selective use of classes of Hagedorn wave packets with specific 
geometries could prove useful for designing meshfree numerical discretizations 
of evolution equations with a priori known symmetries. 
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Figure 2: Intensity plot of the absolute value of the two examplary two-dimensional Hagedorn 
wave packets (p® g[Zi](x) (left) and <P 7 'q[Z 2 ] (rr) (right), e = Hr 1 . Darker colouring represents 
higher absolute values. 



Figure 3: Intensity plot of the absolute value of the two-dimensional Hagedorn wave packets 
ifg c\Zz\(x) (left) and <p§ 7 [Z 2 , Z$](x) (right), where e = Id -1 . 
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Appendix A. Proof of Proposition [4] 


Using the found structure of the polynomials and the corresponding ladder 
operators, we can give a simple proof of Proposition [4] 


Proof. We prove the assertion by induction over \k\. We know that the raising 
operator A I from <© only creates polynomial prefactors in front of The idea 
is to rewrite A I as an operator acting on the prefactors only, and identifying it 
with the raising operator of the polynomials. 

Since the case k = 0 is trivial, we assume the assertion to be true for some 
k £ Then, for j = 1,..., d, by © we have 


ip% +e .[Z,Y\(x) 


y/(k + ej)\ 


= - 7 =L=A]\Y}'p%[Z,Y\(x) 

V k J + 1 

[X*(ieV x -PQ- 1 x)+K*x].q}f ip%[Z\{x). 


Hence, the result is true as long as 

Qk+e, = [—y/eX*W x + ^B*Q~'x\ jq f (^B*Q~ l x) , 


which follows by invoking the polynomial raising operator of Lemma [6j It is 
easy to see that M is symmetric since the components of A^ commute by the 
isotropy condition. 

The structure of the matrix M follows from the fact that both Q~ X Q and 
QQ* are symmetric and 


B*Q* = % -(K*Q - X*P)Q* = 1 -{K*QQ - X*PQ T - 2 iX*) 
since PQ* = (QP* + 2iId) T . Hence, 


B*Q*Q~ t B = l -(K*QB - X*PB) + X*Q~ t B = l -Y*PLZB + X*Q~ T B 

= ^Y*nzz T nY + x*q~ t b = -^y*GzJ" + x*q~ t b, 


where the last equality is due to the isotropy of Y. 


□ 
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Appendix B. Tensor product representation 

By recursively applying Proposition [8] one can derive an expansion of the 
general polynomials qjf 1 in terms of tensor products of univariate Hermite poly¬ 
nomials associated with the diagonal entries of M. Moreover, the required 
number of summands depends only on the number of offdiagonal entries of M 
for which M, 3 ^ 0, and the corresponding indices fcj, and kj. 

Proposition 13 (Tensor product representation). Let M £ C dxd be symmetric, 
and suppose that there are exactly n < d{d — l)/2 different off-diagonal index 
pairs 1 < ay < /3j < d, j = 1 for which M a ,p. = A 3 0. Then, for 

k£N d , 

fc"M= E (b.i) 

^ gN n V /V / i= l 

e j <min{k aj ,k f 3 j } 

with standard multiindex notation, e.g. k a £ N n with ( k a ) 3 = k aj . The index 
matrix E £ N dxn is defined by 


Eij (c ot, T )j ■ 


(B.2) 


Proof. We start by recalling (121, which can be rewritten as 


Qk 


min {kn^kp.} 

0 ) = J2 m] 


m—0 


m J \ m 


(B.3) 


for all j = 1 ,,n. One can use the matrix E in order to write 


k - m{e aj + ep } ) = k - ( Eme 3 ) 

where ej denotes the j-th unit vector in R n . Iterating this procedure until all 
offdiagonal entries of M are deleted, completes the proof. □ 
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